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ABSTRACT 

How do T Tauri disks accrete? The magneto-rotational instability (MRI) supplies one means, but 
protoplanetary disk gas is typically too poorly ionized to be magnetically active. Here we show that 
the MRI can, in fact, explain observed accretion rates for the sub-class of T Tauri disks known as 
transitional systems. Transitional disks are swept clean of dust inside rim radii of '--^lO AU. Stellar 
coronal X-rays ionize material in the disk rim, activating the MRI there. Gas flows from the rim to the 
star, at a rate limited by the depth to which X-rays ionize the rim wall. The wider the rim, the larger 
the surface area that the rim wall exposes to X-rays, and the greater the accretion rate. Interior to the 
rim, the MRI continues to transport gas; the MRI is sustained even at the disk midplane by super- keV 
X-rays that Compton scatter down from the disk surface. Accretion is therefore steady inside the rim. 
Blown out by radiation pressure, dust largely fails to accrete with gas. Contrary to what is usually 
assumed, ambipolar diffusion, not Ohmic dissipation, limits how much gas is MRI-active. We infer 
values for the transport parameter a on the order of a percent for the prototypical systems GM Aur, 
TW Hyd, and DM Tau. Because the MRI can only afflict a finite radial column of gas at the rim, disk 
properties inside the rim are insensitive to those outside. Thus our picture provides one robust setting 
for planet-disk interaction: a protoplanet interior to the rim will interact with gas whose density, 
temperature, and transport properties are definite and decoupled from uncertain initial conditions. 
Our study also supplies half the answer to how disks dissipate: the inner disk drains from the inside 
out by the MRI, while the outer disk photoevaporates by stellar ultraviolet radiation. 

Subject headings: accretion, accretion disks — X-rays: stars — stars: pre-main-sequence — solar system: 
formation — magnetohydrodynamics — planets and satellites: formation 



1. INTRODUCTION: TRANSITIONAL DISKS AND THE MRI 

■rj" . Transitional disks surrounding pre-main-sequence, 
C^l solar-type stars are identified by their spectral energy 
distributions (SEDs). At 1-10 microns wavelength, their 
\^ , excesses are weaker than those of classical, flat-spectrum 
■ T Tauri disks, while at longer wavelengths, their fluxes 
f — ' are as strong as those of any classical T Tauri SED. The 
deficit at infrared wavelengths is consistent with tran- 
;V sitional disks having large inner "holes" that are fairly 
. \ transparent in the continuum. Figure 1 illustrates the 
. situation. Outside the hole radius aiim, optically thick 
^ dust abounds. Inside arim, only trace amounts of opti- 
ca \ cally thin dust are present. 

Some of the more well-known transitional systems are 
listed in Table 1; these include GM Aur (arim ~ 24 AU; 
Calvet et al. 2005, hereafter COS), TW Hyd (a^m ~ 1-4 
AU; Calvet et al. 2002, hereafter C02; Hughes et al. 2007; 
Ratzka et al. 2007), DM Tau (anm ~ 3 AU; COS), and 
CoKu Tau/4 (anm « 10 AU; D'Alessio et al. 2005). In- 
terferometric measurements verify the existence of in- 
ner disk holes (Hughes et al. 2007; Ratzka et al. 2007; 
J. Brown, personal communication). Transitional disks 
may bridge the evolutionary gap between conventional 
T Tauri disks that do not have large inner clearings and 
debris disks that are entirely optically thin. Rim radii 
are too large to represent dust sublimation fronts, fu- 
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eling speculation that inside the hole, either (a) grains 
have grown too large for the medium to remain optically 
thick, (b) planets have consumed or torqued away ma- 
terial, and/or (c) the inner disk has drained onto the 
central star by accretion. This paper concerns possibil- 
ity (c) — or more accurately, how accretion in transitional 
disks is ongoing. 

Transitional disk holes are not empty. Though nearly 
devoid of dust, the region interior to Orim often contains 
accreting gas. Near-ultraviolet excesses, when observed, 
imply stellar accretion rates of M ^ lO~^M0yr~^ 
(MuzeroUe et al. 2000; C02; Bergin et al. 2004; COS); see 
Table 1. We seek here to explain the origin of such disk 
accretion. We propose that gas is leached from the inner 
disk rim by the magneto-rotational instability (MRI), a 
linear instability that amplifies magnetic fields in out- 
wardly shearing disks and drives turbulence (Balbus & 
Hawley 1998). X-rays emitted by the hyperactive corona 
of the central star irradiate and ionize the disk rim, ac- 
tivating the MRI there. Magnetically active gas diffuses 
inward from the rim to the star. The inner hole grows as 
the MRI eats its way out. The accretion rate interior to 
the ever-expanding rim is steady (constant with radius) 
and entirely set by the accretion rate established at the 
rim. 

Historically, whether the MRI is relevant to T Tauri 
disks has been a prospect fraught with doubt (see, e.g., 
Hartmann et al. 2006), as protoplanetary disk gas is typi- 
cally too cold and/or too dusty to couple well to magnetic 
fields. The following two conditions must be met for the 
MRI to be viable. First, magnetic fields must be frozen 
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Fig. 1. — Schematic of a transitional protoplanetary disk accret- 
ing by the MRI. Transitional disks have large rims (arim ~ 10 AU) 
inside of which dust is largely absent. Hard X-ray radiation from 
the young stellar corona photoionizes dusty rim gas and heats it 
to a temperature T* . Only a limited gas column A'^* is sufficiently 
ionized to be MRI-unstable and drain inwards. Stellar radiation 
pressure clears infalling gas of dust. 

into whatever plasma is present. To freeze magnetic flux, 
the magnetic Reynolds number must be sufficiently large: 

^^^^-^(au) (took) (io^)>^^- 

(1) 

Here we follow Fleming et al. (2000) in using the gas 
sound speed Cg and the vertical density scale height 
h = Cg/ft for characteristic velocity and length scales, 
respectively; il is the Kepler orbital frequency; D is the 
magnetic diffusivity; T is the gas temperature; a is the 
disk radius; and a^c is the fractional electron density. 
Numerical simulations by Fleming et al. (2000) suggest 
Re* ~ 10^-10^, depending on the geometry of the seed 
field (but see Pessah, Chan, & Psaltis 2007 for a critical 
re-analysis of results related to the zero-net flux case; it 
is possible only the net vertical and net toroidal cases 
have their values of Re* ~ 200-3000 well measured) . 

A second criterion — more often than not neglected in 
the protoplanetary disk literature — requires that neutral 
hydrogen molecules, constituting the bulk of disk mat- 
ter, be dragged inward by accreting plasma (Blaes & Bal- 
bus 1994; Mac Low et al. 1995; Hawley & Stone 1998; 
Kunz & Balbus 2004). That is, ambipolar diffusion of 
the magnetized plasma out of the (overwhelmingly) neu- 
tral medium must not defeat the MRI. To couple the 
neutrals to the ions, a given H2 molecule needs to collide 
with enough ions within ~l/r2, the e-folding time of the 
instability: 

Am = > Am* , (2) 

where (3in ~ 1-9 x 10~^ cm'^ s^^ is the coUisional rate co- 
efficient for ions to share their momentum with neutrals 
(Draine, Roberge, & Dalgarno 1983) and n is the number 
density of hydrogen molecules (so that Xcn is the num- 
ber density of ions). Linear growth rates for the MRI 



drop dramatically with the degree to which criterion ([2]) 
is not satisfied (Kunz & Balbus 2004) . The critical num- 
ber Am* of neutral-ion collisions is measured in numer- 
ical simulations to be about 10^ (Mac Low et al. 1995; 
Hawley & Stone 1998). To our knowledge, these pioneer- 
ing simulations have not been followed up or verified at 
higher resolution; therefore we regard this measurement 
of Am* , like the measurement of Re* , with caution. Cri- 
teria llj and 12]) are more easily met the more disk gas 
is ionized; criterion ([T]) depends on the fractional ioniza- 
tion, whereas ([2]) depends on the total ion density. 

Transitional disks furnish an excellent first laboratory 
in which to investigate whether the MRI plays a role, if 
any, in how T Tauri disks accrete. Transitional disk holes 
are less plagued by sub-micron-sized dust that threatens 
to scour the region of plasma. Furthermore, the rim wall 
provides an especially simple geometry for analyzing how 
X-rays ionize disk material: if the wall is not severely 
shadowed from the star by material interior to the rim 
(see ij5.2p . one need only treat the 1-D radial transfer of 
radiation perpendicular to the rim wall. 

X-ray driven MRI is not a new idea — it was originally 
introduced by Glassgold, Najita, & Igea (1997, here- 
after GNI97) — but our study applies it for the first time 
to transitional disks and to identify the correct crite- 
rion for where the MRI can operate under these condi- 
tions. We argue that the entire region interior to the 
rim, including the midplane, is MRI-active. This con- 
trasts with previous studies of X-ray driven (GNI97; Igea 
& Glassgold 1999, hereafter IG99) or cosmic-ray driven 
(Gammie 1996, hereafter G96; Sano et al. 2000) MRI 
in which accretion is confined to the surface layers of 
a disk whose midplane properties cannot be calculated 
from first principles. Moreover, such surface accretion is 
unsteady (G96). The accretion rates that we derive here 
are steady and can be directly compared to observations. 
Thus, transitional disks not only provide insight into the 
late stages of disk evolution, but also offer the first clean 
application of ideas — namely, MRI triggered by non- 
thermal ionization processes (G96; GNI97) — pioneered 
for more complex systems. 

In [J2I we construct a simple model for the rate of 
mass flow from the X-ray illuminated rim. We solve 
the equations of ionization and energy balance to decide 
how much of the rim satisfies criteria H]) and ^ and 
accretes. In addition to accounting for the usual disso- 
ciative recombination of molecular ions with free elec- 
trons (e.g., GNI97), our treatment of ionization balance 
includes recombination onto dust grains (Umebayashi & 
Nakano 1980) and charge transfer to free atomic metals 
(Oppcnheimcr & Dalgarno 1974, hereafter OD74; Fro- 
mang, Terquem, & Balbus 2002). We show how criterion 
([2]), neglected in many studies of protoplanetary disks 
(e.g., G96; GNI97; IG99; Sano et al. 2000; Fromang et 
al. 2002; Matsumura & Pudritz 2003), supersedes crite- 
rion (U)). We establish a relation between accretion rate 
M and rim radius Orim and compare to observations. 

In we demonstrate how the MRI continues to trans- 
port gas inward at radii a ^ Orim • Though the midplane 
of the interior disk is optically thick to incident X-rays, 
Compton scattering permits enough X-rays to reach the 
midplane (IG99) that the MRI operates there as well. 
Accretion is therefore steady from the rim to the star. 
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In 21 we study in a preliminary way how stellar radia- 
tion pressure blows dust grains out of accreting gas and 
can keep the disk inside the rim relatively transparent. 
We make a connection to the small but non-zero opti- 
cal depth inferred for the inner disk from SED modelling 
(C02; COS) to constrain the abundance of micron-sized 
grains in the rim. 

Our work leaves many issues unresolved and calls at- 
tention to outstanding questions associated with the 
MRI. These are discussed in after a summary of our 
principal findings. 

2. ACCRETION FROM THE RIM 

Ionization of rim gas is maintained by stellar X-rays.'^ 
As measured by Chandra, the X-ray spectra of T Tauri 
stars in Orion often have separate soft and hard compo- 
nents, modelled as emission from thermal plasmas hav- 
ing characteristic energies kT w 1 keV and kT w 3 keV, 
where k is Boltzmann's constant (Wolk et al. 2005, see 
their Figure 12). The luminosities of the two components 
are comparable, and depending on the star, range from 
10^^ to lO^^ergs"^. Soft, sub-keV X-rays are thought 
to be emitted by shock-heated gas in accretion columns 
(Kastner et al. 2002; Stelzer & Schmitt 2004), while hard, 
super-keV photons may arise from magnetic flares in stel- 
lar coronae (e.g., Wolk et al. 2005). 

As the penetrating power of X-rays increases rapidly 
with photon energy, we are primarily interested in the 
hard, presumably coronal emission. Unfortunately, we 
are unable to locate detailed X-ray spectra at photon en- 
ergies Ex ^ 3keV for the prototypical transitional sys- 
tems discussed in this paper (GM Aur, TW Hyd, DM 
Tau, and CoKu Tau/4). What we glean from the litera- 
ture is summarized in Table 1 . In the absence of detailed 
spectra, we assume for our computations below that each 
source emits Lx — lO^^ergs"^ at a characteristic en- 
ergy Ex = 3 keV. This assumption is supported circum- 
stantially by the aforementioned Chandra observations of 
hard X-ray emission from Orion prc-main-sequence stars. 
It is, furthermore, consistent with our findings in Table 
1, even for TW Hyd, whose spectrum is known to be 
atypically soft (Stelzer & Schmitt 2004). 

Stellar X-rays penetrate the exposed rim wall to a hy- 
drogen column density N (see Figure 1). A limited col- 
umn N* will satisfy criteria (|T|) and ([2]). The MRI-active 
rim contains mass Mrim ~ AT:N*arimhfi, where h — Cs/f2 
is the vertical density scale height, Cs = {kT* / /i)^/^ is the 
gas sound speed, is the Kepler orbital frequency, T* 
is the kinetic temperature of MRI-active gas at the rim, 
and /X « 3 X 10~^^g is the gas mean molecular weight. 
This mass flows from arim to '^arii„/2 over the diffusion 



time ^diff' ^ 
diffusivity: 



^/v, where v = acgh is the turbulent 



M 



3Mn 



127raiV*a2,„(fcr*) 



^3/2 



(3) 



^ Galactic cosmic rays at GeV energies yield ionization rates of 
at most ~10~^'^s~^ per H2 (Spitzer 1978, p. 116). Actual cosmic 
ray ionization rates in T Tauri disks are likely much lower because 
of shielding by magnetized T Tauri winds. Even the much weaker 
solar wind modulates the cosmic ray flux at Earth by as much 
as ^10% with solar cycle (Marsh & Svensmark 2003). We ignore 
cosmic ray ionization; it is easy to check that their contribution is 
negligible compared to ionization by stellar X-rays. 



where G is the gravitational constant and the extra fac- 
tor of 3 follows from a more accurate derivation (using 
M ~ SnT,!/, where S is the disk surface density; Frank, 
King, & Raine 1992). As measured by numerical sim- 
ulations of the MRI, the transport parameter a might 
be ~0.01-0.3, depending on the strength and geometry 
of the background field (Fleming et al. 2000; Sano et 
al. 2004; Pessah et al. 2007; see also ^5?2)) . Notice how 
M in equation ([3]) does not depend explicitly on the sur- 
face density profile S(a) of the disk. Given observations 
of Orim and M*, one need only compute the rim-specific 
variables TV* and T*, which we do below. 

2.1. Ionization Balance 

Determining N* requires that we calculate the degree 
of ionization as a function of the radial column TV pen- 
etrated by X-rays at the rim. Per H2 molecule, the ion- 
ization rate is 



c 



riaxexp{-Nax) 



(4) 



where ctx » 4 x lO-^^{Ex/SkeV)-^-^'^ cm'^ is the pho- 
toionization cross-section per II2 (assuming solar abun- 
dances) and r/ « 81 accounts for the number of sec- 
ondary ionizations produced per absorbed 3-keV photon 
(GNI97). 

Freshly ionized rapidly converts to molecular ions 
such as HCO"*". Most of these molecules dissociate in col- 
lisions with electrons, while some transfer their charge 
to gas-phase atomic metals such as magnesium (Fro- 
mang et al. 2002; OD74). The fractional number abun- 
dance a;inoi+ of molecular ions relative to neutral hydro- 
gen molecules is given by the equilibrium relation 



(5) 



where Pdiss ~ 3 x 10-^(r/230 K)-i/2 cm^ s^^ and A w 
10~® cm'^ s~^ are rate coefficients for dissociation (Glass- 
gold, Lucas, & Omont 1986) and charge transfer (OD74), 
and Xo and a;niot are the fractional number densities of 
electrons and free neutral metals, respectively. To re- 
late the number density n of neutral H2 molecules to N, 
we spread all 47r7Varim'i molecules comprising the X-ray 
penetrated rim into the volume 2Tra'^^^^^h interior to the 
rim to estimate that n ~ 2N / a-^ira- 

Charge transfer to free metals is important for sustain- 
ing the MRI because ionized metals tend to keep their 
charge, neutralizing slowly either by radiative recombina- 
tion with electrons or collisions with negatively charged 
dust grains (OD74; Umebayashi & Nakano 1980; Glass- 
gold et al. 1986): 



(6) 



where a;„iet+ is the fractional abundance of metal ions, 
^rec ~ 4 x 10-i2(T/230K)~i/2cm3s^i is the radiative 
recombination coefficient, and f3gr is the recombination 
coefficient, measured per II2, for grains. We assume the 
total grain surface area available for recombination is 
dominated by grains of radius s containing a fraction 
Zs of the total mass in gas and dust. It follows that 
/3gi. w 3 X 10-20(^m/s)(Z,/10-'*)(r/230K)i/2cm3s-i. 
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TABLE 1 

Properties of Transitional Disk Systems 



Star 






M 


M 


r a 

-'^X.obs 


X-ray band 


X-ray 




(AU) 


refs. 


(lO-'JAfoyr-i) 


refs. 


(ergs-1) 


observed (keV) 


rcfs. 


GM Aur 


24 


1 


2.7-10 


1,2 


1030 


0.2-2.4 


3 


TW Hyd 


1-4'' 


4,5,6 


0.3-0.7'= 


4,7 


1.4 X 10^" 


0.45-2.25 


8 


DM Tau 


3 


1 


0.8-2 


1,2 


2 X 1029 


0.3-10 


9 


CoKu Tau/4 


10 


10 


< 0.1 


11 


NA 


NA 


12 



References. — 1. COS; 2. Bergin et al. 2004; 3. Strom et al. 1990; 4. C02; 5. Hughes et 
al. 2007; 6. Ratzka et al. 2007; 7. MuzeroUe et al. 2000; 8. Stelzer & Schmitt 2004; 9. Guedel 
et al. 2007; 10. D'Alessio et al. 2005; 11. Najita et al. 2007; 12. Konig et al. 2001 

In modelling the hard X-ray spectrum incident upon the rim, we do not simply use the 
observed X-ray luminosity Lx obs, which is typically measured for soft (Ey^ < 3keV) photons 
having low penetrating power. Instead, we assume in our calculations of ionization and thermal 
balance that each star emits Lx = lO^^ ergs~^ at a characteristic photon energy i?x = 3keV. 
^ Controversy exists whether Orim ^ 4AU (C02; Hughes et al. 2007) or is nearly ~1 AU 
(Ratzka et al. 2007) for TW Hyd. For our analysis we adopt the former result, since it is based 
on methods similar to those applied for GM Aur, DM Tau, and CoKu Tau/4. The conclusions 
of our paper are unaffected by this choice. 

Taken to span a factor of 2 to reflect U-band variability (MuzeroUe et al. 2000). 



This estimate, which assumes that every grain is nega- 
tively charged, agrees with that of Umebayashi & Nakano 
(1980), who assume s « 0.1 /.tm and Zg k, 10~^. Though 
unknown, the factor Zg is likely to be considerably less 
than the maximal value permitted by solar abundance 
gas, Zq « 10~^, because of grain growth.'' We take as 
standard values Zg = 10^^ and s = 1 /im. A constraint 
on Zg/ s from observations is contained in ^jJlon radiation 
blow-out. 

In assuming a single grain size, we neglect effects intro- 
duced by smaller-sized grains. The surface area for ion 
recombination may be dominated by very small grains. 
Furthermore, because grains efficiently adsorb electrons, 
they can, if sufficiently numerous, completely clear the 
gas of free electrons and catastrophically reduce the elec- 
trical conductivity. Indeed, this is the case in the Earth's 
lower atmosphere, where charges are carried by aerosols. 
Fortunately, as we describe below in ^2.'6[ our results 
enjoy some margin of safety from these effects. The frac- 
tional electron densities that we derive, 10~^-10~^ (com- 
pared, say, with terrestrial ionization fractions ~10~^^ at 
sea level) , are large enough that the chemistry is arguably 
driven by electrons (but see ji5.2p . 

Charge and number conservation read, respectively. 



''mol+ 



-|- X 



met+ 



and 



'^met,tot — '^met+ ~r ^mct ■ 



(7) 



(8) 



Because free metals are depleted onto grains, Xmct.tot is 
realistically at most ~3 x 10~^, corresponding to 1 out 
of every 10 metals in the free atomic phase (OD74; Ume- 
bayashi & Nakano 1980). We consider values for a;,„et,tot 
ranging from zero (complete depletion) up to 10^^. 

We combine equations ([5]), ([H]), (O, and ([5]) to derive, 
under certain approximations, a quartic equation for 
in terms of C,/n, Xmct.tot, and various rate coefficients. 

* If grain growth were so extreme that an individual grain 
shielded itself from incident X-rays, then we would have to reduce 
the X-ray cross-section crx (which is measured per H2). 



Equations 



|), and jll) yield 



c c 



0. (9) 



n n \Xcf3icc + /3g 
The first term on the left-hand side dominates the second 
term when 

1/2 



A. 



"^met.tot 



'diss 



9 /W^x / _T_ 
V 10-6 J \230K 



(10) 

This is always the case for MRI-active rim material, as 
will be evident below (see Figures 2 and 3). Therefore 
we drop the second term in ([9]) to write 



n I 2 '^/^diss 



1 (/3rcc2^c + /3gr) 



(11) 



Combining equations (O, ([51), and ([5]), we have 

(/^diss*^c ~^ /^t-^mct )'^mct,tot 

which simplifies, according to the same approximation 
embodied in pU)). to 



/^diss'^mct,tot'^e — *^mct 



/^diss'^c 



(13) 



n \(3rccXc + Pi 
We solve (|13|) for a;inct and substitute into (fTTj) to write 

/5diss/^rcc*^c ~^ /^diss/^gr'^o ~^ n Pi'cc') -^q 

-^(/3gr+/3tW,tot)Xe- (^)'^ -0. (14) 

It is safe to ignore Prcc in comparison to Pt , so that 

~ {nL.sil) ( + ^mct.tot) (15) 

- 

This final equation is a quartic and not a cubic (OD74) 
because we allow for the possibility that nearly all gas- 
phase atomic metals might be ionized, via equation ([8]). 
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2.2. Thermal Balance 

Gas is heated by fast photoelectrons at a rate per H2 
molecule of 



(16) 



where / ~ 0.5 is the fraction of X-ray energy deposited 
into heat (Glassgold, Najita, & Igea 2004). Gas cools by 
ro- vibrational transitions of CO at a rate per H2 molecule 
of 



A =3.3 X IQ- 



-25(rp/r^\l/2 cxp(-4292K/T) 

[-^ / l-cxp(-3084K/T) ^ 



a;co("-/cm 



erg s 



(17) 



(Glassgold et al. 2004, their equations C7-C8). Here 
a^co ~ lO^'' is the fractional abundance of CO, whose 
level populations are assumed thermal. Other cool- 
ing channels — dust-gas collisions and Lyman a emission 
from hydrogen — are less important and ignored. We re- 
place n in ((TT)) with our estimate 2iV/arim (see i|2.ip . set 
r = A in thermal balance, and solve for T as a function 
of N. This result is necessary for solving the quartic (jisp 
for Xc because rate coefficients depend on T. 

For our standard parameters for CM Aur (Lx — 
10^9 erg s-\ Ex = 3keV, anm = 24 AU), we find that 
T drops from about 340 to 200 K as increases from 
10^^ to 10^^ cm^^. Results for T are very similar for TW 
Hyd and DM Tau. The insensitivity to input parameters 
is due to the exponential sensitivity of the cooling rate 
A to T. 

2.3. The Active Column and Resultant Accretion Rate 

Taking the run of T with N from our simple model 
of thermal balance in ^2.2\ we solve the quartic (fT5|) nu- 
merically. The solution Xc is displayed against N in Fig- 
ure 2 for parameters appropriate to CM Aur, for various 
choices of Xmet,tot- Evidently the presence of free atomic 
metals raises the ionization fraction by up to two orders 
of magnitude. Results for TW Hyd and DM Tau are 
similar to within factors of two. 

From Xc{N) we compute Am(A^), as shown in Fig- 
ure 3 for the case of CM Aur. If Am > Am* = 100 
strictly determined the viability of the MRI, then we 
would equate N* with the maximum value of N for which 
Am = 100. However, given uncertainties in the value of 
Am* and in the sharpness of criterion ^ in predict- 
ing the non-linear, saturated outcome of the MRI, only 
an order-of-magnitude estimate for TV* 5 x lO^'^cm^^ 
(the value for which Am peaks) seems justifiable. We 
adopt this same value of N* for TW Hyd and DM Tau, 
for which very similar results for Am(A'^) obtain. For 
N — N* = 5 X lO^'^cm"^, we find from our thermal 
balance model that T = T* 230 K. 

The ionization fractions, Xc ^ 10^*-10^^, correspond- 
ing to ~ N* are so large that magnetic Reynolds 
numbers. Re ~ 10^-10*, far exceed the critical values. 
Re* ^ 10^-10^, seemingly required for magnetic flux 
freezing (Fleming et al. 2000). Therefore accretion by 
the MRI in transitional disks is limited by ambipolar 
diffusion, not by Ohmic dissipation. 

Electron densities exceed grain densities by 8-9 orders 
of magnitude, for our standard Zg = 10^** and s = 1 /im. 
Our neglect of a population of aerosols which can scour 
the gas of electrons is arguably safe. The differential 



-1 1 I I I I I l| 1 1 I I I I llj 1 TT 

GM Aur 

a,i„ = 24 AU : 




N (cm"^) 



Fig. 2. — Ionization fraction Xc against column density A'^ pen- 
etrated by X-rays at the rim, evaluated for GM Aur (assuming 
Lx = 10^^ ergs-l and Ex = 3keV). Curves are computed using 
the quartic equation JTSj for x^. Results for TW Hyd and DM Tau 
are similar. Solid, long-dashed, dashed, and dotted curves corre- 
spond to various choices for the gas-phase atomic metal abundance 
a;mot,tot: 10~^, 10"'', 10~*, and 0, respectively. An abundance 
a;mct,tot = 10~^ corresponds to 1 out of 30 metals in the atomic 
gas phase. Free atomic ions, which recombine slowly, can increase 
the ionization fraction by up to two orders of magnitude. 

grain size distribution, if it scaled as s~^-^, would have 
to extend from s = 1 fim all the way down to several 
angstroms for the total grain density to be comparable 
to the electron density. Our findings are also insensitive 
to the uncertain grain recombination coefficient /3gr, be- 
cause metallic ions recombine more readily with electrons 
than with grains. For example, increasing /Jgi- by a factor 
of 30 above our nominal value (such an increase would 
correspond to the same grain size distribution mentioned 
above) decreases the peak value of Am in Figure 3 by a 
factor of 2, from 140 to 70. 
Armed with N* and T* , we use equation ^ to plot 

M against arim in Figure 4, adjusting a as necessary 
to reproduce observed accretion rates for GM Aur, TW 
Hyd, and DM Tau. The offsets between the different lines 
drawn for different systems arise mostly from differences 
in M* and a, since N* = 5 x 10^"^ cm~^ is held fixed and 
r* (TV*, arim) varies little. The slopes of the lines (M cx 
"^rim' approximately) reflect essentially the increase in the 
rim wall's surface area with Orim- 

Best-fit a's are 0.007-0.035. Such values are of order 
those seen in numerical simulations of the MRI, depend- 
ing on the strength and geometry of the background field 
(Fleming et al. 2000; Sano et al. 2004; Pessah et al. 2007; 
see also 



^ Intriguingly, our best-fit o's resemble those required to fit the 
light curves of dwarf novae during their "cold" states (a 0.02— 
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Fig. 3. — Ambipolar number Am against column density N 
penetrated by X-rays at the rim, evaluated for GM Aur (assuming 
Lx = lO^^ergs-l and Ex = 3keV). Results for TW Hyd and 
DM Tau are similar. Solid, long-dashed, dashed, and dotted curves 
correspond to various choices for the free atomic metal abundance 
J^metjtot: 10~®, 10~^, 10^*, and 0, respectively. An abundance 
aimet.tot = 10~^ corresponds to 1 out of 30 metals in the atomic 
gas phase (i.e., 29 out of 30 depleted onto grains). The ambipolar 
numbers required for accreting plasma to entrain neutral hydrogen 
are measured in numerical simulations to be Am* ~ 10^ . Roughly 
speaking. Am ~ Am* for A'" = N* ~ 5 X 10^^ cm~'^, provided 

The transitional system CoKu Tau/4 is not detected 
in X-rays (Konig et al. 2001) and has an unmeasur- 
ably small accretion rate (M < 10^^'-' Mq yr~^; Najita et 
al. 2007; D'Alessio et al. 2005), facts that at face value 
are consistent with our theory. Possibly CoKu Tau/4 
has fewer free atomic metals or — and this is more readily 
testable — a softer X-ray spectrum than the other sources. 
If, say, Am < 10 at its rim, then we estimate its 3-keV 
luminosity Lx < 102^{10-^/a;inet,tot) ergs'^ 

We expect little accreting material outside aj-im- The 
X-ray heated rim has a vertical thickness greater than 
that of material immediately outside. Thus, as depicted 
in Figure 1, some fraction of the disk beyond the rim 
will dwell in the rim's X-ray shadow and be magnetically 
inert. Whatever distant material lies outside the rim's 
shadow and might be accreting from X-ray irradiation 
is ignored. See also ^ 35. II for a reason to doubt that the 
outermost regions of disks accrete. 

3. STEADY ACCRETION IN THE INTERIOR DISK 

0.04; Lasota 2001). This resemblance seems sensible insofar as 
white dwarf accretion disks during their quiescent phases and pro- 
toplanetary disks both contain relatively cold, poorly ionized gas. 
However, Gammie & Menou (1998) argue that dwarf nova disks in 
quiescence do not support the MRI. These authors consider the rel- 
evant magnetic Reynolds numbers. Re ~ 3700, too low. Our own 
view is that such Reynolds numbers do not proscribe a weakened 
form of the MRI (see also Fleming et al. 2000). 



Fig. 4. — Accretion rates M versus rim radii a^im. Points with 
error bars represent observed data for GM Aur, TW Hyd, and 
DM Tau (see Table 1 for references). Solid lines are calculated 
according to equation l(3)l, for N* = 5 X 10^^ cm"'^, Mt, as given 
by the literature for each star, and T* as computed as a function 
of a^ini and Af« for fixed N*. Offsets between lines mostly reflect 
differences in M* and a. The accretion rate M increases with arim; 
the larger the hole, the more surface area the rim wall exposes to 
X-rays. The transport parameter a labelling each curve is chosen 
to reproduce the observations. Fitted values of a range from 0.007 
to 0.035, close to values reported in simulations of the MRI. 

Once dislodged from the rim, gas must still travel up to 
three decades in radius, from ~10 to ~0.01 AU, to reach 
the stellar surface. We now argue that X-ray-driven MRI 
continues to provide transport at a ^ Orim- At ~0.01 
AU, the MRI should be well sustained by thermal ion- 
ization. It remains to decide whether criteria ([T|) and ^ 
are satisfied at the midplane at radii of ~0.1 to ~1 AU. 
We need to estimate the temperature and surface density 
of material there. 

Observed spectra demand that inside the rim, the disk 
contain so few grains as to be optically thin at mid- 
infrared wavelengths (C02; C05). We can understand 
this as a consequence both of the limited number of 
grains contained within the MRI- active mass Mrim, as 
well as stellar radiation pressure, which blows out a large 
fraction of grains having sizes s < 1/zm (see Sj4]). Hav- 
ing lost its primary source of continuum opacity, the gas 
heats not by incident starlight, but by accretion (Ohmic 
dissipation). Cooling proceeds through line emission, at 
a rate that is difficult to estimate because it requires solv- 
ing simultaneously for the thermal, chemical, and exci- 
tation state of the gas. 

We do not attempt to construct a full-blown thermal 
model of electrically resistive, line-emitting gas. Instead, 
we find it adequate to estimate the temperature of mid- 
plane gas at a ^ in terms of its minimum possible 
value, obtained by equating the energy flux from accre- 
tion with that emitted by a blackbody: 
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T« 50 



M 



1/4 



/ a \-3/4 



10-9Moyr-i / VAU/ 



TK. 



(18) 



where T > 1 parameterizes our uncertainty. Mass con- 
tinuity then implies a vertical hydrogen column density 
of 



iV_L « 5 X 10^4 



M 



\ 3/4 



-3/4 



(19) 

to the midplane. We have normalized a to the average 
of values obtained for GM Aur, TW Hyd, and DM Tau 
(Figure 4). 

Encouragingly, equation (fT9|) yields surface densities 
that are entirely consistent with those inferred by Salyk 
et al. (2007) for GM Aur and TW Hyd: at a - 0.2- 
0.5 AU, these authors find N± ~ 0.3-3 x lO^^cm^^ (H2 
surface densities of 1-10 gcm~^ assuming xqo = 2 x 
10"'^), based on observations of rovibrational emission 

from CO. These same observations suggest that T « 3 
(see their Figure 2), and we normalize our results below 
to this value. 

To estimate Re and Am at the midplane, we first ob- 
tain the ionization rate C from IG99, who calculate by 
Monte Carlo methods how X-rays incident upon disk 
surfaces propagate vertically and radially through pro- 
toplanetary disks. Their calculations account for mul- 
tiple Compton scattering off bound electrons and for 
a thermal spectrum of photon energies. According to 
their Figure 3, for (a) a kT — 3keV plasma emitting 
Lx = 10^^ erg s^^, (b) N± as given by equation p^ . 
and (c) T = 3, the ionization rates at the midplane equal 
C ~ 10-16 s-i at a = 1 AU and C ~ 2.5 x lO-^S"! at 
0.1 AU. 

We insert these rates into the quartic equation (fT5|) for 
Xc, setting /3gr = because dust is largely absent in the 
inner disk. The quartic yields Xc ~ 2x 10^^ at 1 AU 
and a;c ~ 1 X 10"^ at 0.1 AU, for Xmct.tot = 10"^. These 
ionization fractions correspond to Re ~ 2 x 10^ at 1 AU 
and 900 at 0.1 AU, values that satisfy criterion al- 
beit marginally. Finally, taking n = N±/h = N±^/cs, 
we find that Am w 90 at 1 AU and Am w 120 at 0.1 
AU. These values for Am, which marginally satisfy cri- 
terion ([2]), are insensitive to variations in T, which pro- 
duce changes in n that nearly cancel changes in Xe- The 
margins of safety for Re and Am would increase by ac- 
counting for thermal ionization. 

We conclude that the MRI plausibly operates every- 
where interior to the rim, even at the disk midplane, 
carrying steadily inward all of the mass drawn from the 
rim. Note how the MRTactive vertical columns at 0.1-1 
AU, N± ~ 10^^ cm"^, are larger than the MRI-active 
radial column at the rim, N* ~ 5 x 10^"^ cm~^. This fol- 
lows from criterion ([2|), which assigns importance to the 
total density XcU, not the fractional density Xc- Total 
densities increase dramatically from the rim to the star. 

4. RADIATION BLOW-OUT OF GRAINS 

Gas drawn from the rim entrains dust grains. As the 
surface density of inspiralling gas increases towards the 



star, so, too, does the surface density of grains. Left 
unchecked, the increasing concentration of grains would 
render the inner disk optically thick at infrared wave- 
lengths, violating observations. Here we explore how 
stellar radiation pressure can purge grains from gas (see 
also Eisner, Chiang, & Hillenbrand 2006). 

Sub-micron-sized grains around solar-type stars feel an 
outward stellar radiation force that just exceeds stellar 
gravity. Then the time for such grains to travel from a 
to ~2a is 

1 fl/Q' 



tblo 



f'stop 



(20) 



where igtop ^ Pps/{csfin) < 1/Q is the time for grains 
of radius s and internal density pp « Igcm""^ to at- 
tain terminal velocity according to the Epstein gas drag 
law, which applies here because the grain size is smaller 
than the collisional mean free path in gas (e.g., Wei- 
denschilling 1977). Radiation blow-out is faster than 
aerodynamic drift (the latter caused by radial pressure 
gradients in gas; e.g., Alexander & Armitage 2007) by 
{^/^)iim 10^, and so we ignore the latter. 

We compare tbiow to tdiff ^ a^/^^, the time for gas to 
diffuse from a to ~a/2. At the rim of the disk of GM 
Aur, the timescales match by coincidence: ibiow.rim ~ 
^diff.rim ~ 1 X 10^ yr. This indicates that about half of the 
dust in Mj-im is expelled, leaving the other half entrained 
with accreting gas and spread between Orim and ~arim/2. 
The geometric optical depth of entrained dust, measured 
perpendicular to the midplane, is 



2N*iiZs h 
PpS a 



0.5 



10- 



(^) 



(21) 



We can try to check our estimate for Trim against ob- 
servations of the mid- infrared SED. Because radiation 
at a wavelength of 10 /im originates from grains hav- 
ing temperatures of ~300 K, observed 10 pm spectra of 
GM Aur constrain the optical depth in grains at ~1 
AU but not near the rim at 24 AU, where starlight- 
heated grains would have colder temperatures of '--^lOO 
K.6 Therefore we cannot directly compare our calculated 
Trim with 10 ^m observations for GM Aur. 

For the case of TW Hyd, however, we can more easily 
make this comparison, because its rim is located at ^4 
AU. Scaling to the parameters of that system, we find 
ibiowaim 5 X 10^ yr and tdiff,rim 2 x 10^ yr, which 
suggests that somewhat more than half of the dust in 
Mj-im is expelled. Reducing our estimate in (|2ip by an 
additional factor of 2, and accounting for the smaller as- 
pect ratio (ft./a)rim, we estimate Trim ~ 0.1 for TW Hyd. 
Observationally, the vertical optical depth of dust inte- 
rior to the rim of TW Hydra's disk is tiq « 0.05 at a 
wavelength of 10 /im (C02). This is essentially the same 
as the geometric optical depth for micron-sized grains 
because of the silicate resonance band at 10 pm wave- 
length. Thus, our crude estimate of r^m 0.1, based on 
assumed values for Zs = 10^* and s = 1 /im, is within 
a factor of two of the observed optical depth. Such ap- 
proximate agreement we consider acceptable, and helps 
to justify our choice earlier for Zs/s in calculating the 
grain recombination coefficient /3gr ( §2.ip . 

® X-ray heated gas and starlight-heated grains generally have 
different temperatures. 
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To satisfy the observation that the disk remain op- 
tically thin at a <C arim, the order- half of rim dust 
that lies between Orim and must fail to pene- 

trate inside a < arim/2. Radiation pressure ensures 
that dust does not continue to inspiral, provided that 
^blow/^ diff oc riT^/^a^/^ decrease with decreasing a. This 
seems likely to obtain because T should drop sharply just 
inside arim/2 once midplane gas becomes too dense to be 
heated effectively by X-rays. Unfortunately, a conclusive 
statement cannot be made without constructing a model 
that smoothly bridges conditions at a ~ aiim with condi- 
tions at a ^ flrim- This would require knowing in detail 
how gas heats and cools at all radii. 

5. CONCLUDING REMARKS 

We summarize our results in ^b.W and discuss topics 
for future research in ^ 35.21 

5.1. Summary 

How do T Tauri disks accrete? We have answered this 
longstanding question for the sub-class of objects known 
as transitional systems, disks having AU-sized or larger 
inner holes cleared of dust. Super-keV X-rays, emitted 
by the young star as a consequence of its enhanced dy- 
namo, irradiate the inner disk rim. A portion of the 
rim is sufficiently ionized that it becomes unstable to 
the magneto-rotational instability (MRI) and diffuses in- 
ward. The MRI, sustained even at the disk midplane by 
hard X-rays, transports gas steadily from the rim all the 
way to the stellar surface. We infer transport parame- 
ters a on the order of a percent; such values are indeed 
yielded by the MRI, according to contemporary numeri- 
cal simulations. 

We feel the following features of our theory deserve 
emphasis. 

1. Accretion rate vs. hole size. All other factors 
(e.g.. X-ray luminosity, stellar mass, free atomic 
metal abundance) being equal, the accretion rate 
M should increase nearly as the square of the 
rim radius Orim, simply from the increased surface 
area of the rim wall that is exposed to stellar X- 
rays. Late-breaking data (C. Espaillat, personal 
communication) for the transitional disks CS Cha 
(M w 1.2 X lO-**M0yr-i, Orim ~ 43 AU) and SZ 
18 (M w l.TxlQ-^Moyr-^ UHm « 14 AU) appear 
to support this trend of M with Orim as plotted in 
our Figure 4. Of course, this trend is expected to 
apply only for transitional disks for which Orim > 1 
AU (see the last issue raised in 

2. Ambipolar dijfusion vs. Ohmic dissipation. Often 
it is assumed that Ohmic dissipation controls where 
the MRI can and cannot operate. This criterion, 
embodied in equation ([T]) through the magnetic 
Reynolds number Re, assigns importance to a^o, the 
ratio of the density of the most mobile charge car- 
riers (electrons) to that of resistive species (neutral 
H2 molecules). Large values of Re oc Xc imply that 
magnetic flux is frozen into plasma, enabling the 
background shear flow to amplify magnetic flelds. 

We have shown instead that ambipolar diffusion, 
not Ohmic dissipation, dictates the boundaries of 



MRI-active zones. Magnetic flux can be perfectly 
frozen into plasma; yet if the neutrals do not cou- 
ple mechanically to the plasma, the neutrals would 
be immune to the MRI. Protoplanetary disk gas 
is overwhelmingly neutral. Thus, as equation ([2]) 
states, a given neutral H2 molecule must collide 
with as many (momentum-bearing) ions as possi- 
ble within a dynamical time for the MRI to afflict 
protoplanetary disk gas. Ambipolar diffusion lim- 
its the accretion flow from the rim, and is only 
marginally defeated at the midplane far inside the 
rim. 

While we have demonstrated that criterion ([2]) out- 
weighs ([T]) only for transitional disks, we suspect 
the same conclusion holds for conventional T Tauri 
disks as well. For example, consider Figure 7 of 
IG99, which seeks to describe the depth of the 
MRI-active layer in disks without holes. At, say, 
a = 1 AU, at the bottom of the supposedly ac- 
tive layer of their model, the number density is 
n ~ N±/h ~ 2 X 10"cm~3. By IG99's equation 
(23), the ionization fraction there is a;e 4 x lO^^"' 
(for their standard a = 1). It follows that Am = 
XcuPin/^ ~ 10^^. This falls four orders of mag- 
nitude short of the critical Am* ^ 10^ seemingly 
demanded by numerical simulations (Mac Low et 
al. 1995; Hawley & Stone 1998). Therefore we dis- 
agree with IG99's determination of the depth of 
the active layer; their layer does not satisfy crite- 
rion 

Because the ambipolar diffusion criterion ^ as- 
signs importance to the total ion density rather 
than the fractional electron density, it is hardest to 
satisfy in the most rarefied regions of disks, either 
at large vertical heights or large radial distance. 
This is just opposite to what has been concluded in 
the literature, that MRI-active zones are restricted 
to the surfaces and the outermost portions of disks. 

3. Robustness. Disk properties inside the rim are 
insensitive to those outside, since the MRI can 
only draw a radial column of N* 5 x 10^^ cm~^ 
from the rim at any time. In other words, un- 
like the usual situation, estimating M does not re- 
quire specifying the detailed surface density profile 
of the disk S(a); we need only know those rim- 
specific variables in equation ([3]): arim, a, A^*, and 
T*. Our picture therefore provides a robust set- 
ting for theories of how planets form and how their 
orbits evolve (e.g., Goldreich & Sari 2003) within 
transitional disks. A protoplanet lying interior to 
the rim will interact with gas whose density, tem- 
perature, and transport properties are definite and 
decoupled from uncertain initial conditions. 

4. How disks dissipate. Our study also supplies part 
of the answer to how disks dissipate. Except- 
ing matter that gets locked into planets, the in- 
ner disk drains from the inside out by the MRI, 
while material beyond the rim photoevaporates by 
stellar ultraviolet radiation (Alexander, Clarke, & 
Pringle 2006ab). Indeed, the models of Alexan- 
der et al. (2006ab) require the disk inside a few 
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Alls to drain onto the star before ultraviolet radi- 
ation can photoevaporate the outer disk efficiently. 
Their analysis did not specify the means by which 
the inner disk evacuates. Our work identifies the 
mechanism: it is the MRI, leaching material from 
the disk rim and causing the inner hole to grow. 
How long the MRI takes to eat its way out depends 
on the unknown surface density S of the original 
disk. The clearing time tciear ~ ^Orim/-^ ^ 10^ yr 
if S ~ lO^gcm"^ at onm ~ 10 AU, comparable to 
that of the minimum-mass solar nebula. 

5.2. Unresolved Issues 
Numerous areas exist for improvement. 

1. Radiation blow-out. We have established that a 
fraction on the order of 50% of the sub- micron-sized 
grains in the MRI-active rim gets blown back into 
the rim by stellar radiation pressure. The remain- 
ing fraction, which gets entrained to r^a,-i^/2 by 
the accretion flow, ultimately must also be blown 
back to ensure that transitional disk holes stay as 
transparent as observed. To prove this last point 
requires understanding in detail how gas heats and 
cools at all radii from Orim to the stellar surface, 
and checking that ibiow/^diff « uT^/'^a^^'^ decreases 
with decreasing a. 

2. Self- shadowing. Throughout this work we have as- 
sumed for simplicity that the accretion disk interior 
to the rim does not obstruct the rim wall from stel- 
lar X-rays. Yet there must be some self-shadowing, 
which if sufficiently severe will convert the problem 
we have solved to one involving surface irradiation 
and the two-dimensional radiative transfer of X- 
rays (IG99). Even if this were the case, however, 
midplane gas at the rim may still be MRTactive by 
virtue of X-rays that Compton scatter downward 
from the disk surface (IG99; S|3|). 

The rim might be tall enough that a good fraction 
of its lowermost vertical scale height escapes being 
shadowed. In our picture, rim gas is optically thin 
to X-rays, and is therefore well-heated by them (see 
also Glassgold et al. 2004). This source of heating 
is not enjoyed by midplane gas in the inner disk, 
which is optically thick to X-rays. Thus we envi- 
sion a sharp drop in temperature and a concomi- 
tant drop in disk thickness just inside the rim. As 
was the case for radiation blow-out, addressing the 
problem of self-shadowing requires a global treat- 
ment of thermal balance. 

3. The value of a. We have estimated that a's on 
the order of a percent are required to explain the 
accretion rates of GM Aur, TW Hyd, and DM Tau. 
Such values are observed in numerical simulations 
of the MRI under certain assumptions about the 
magnitude and geometry of the seed magnetic field. 
Are these assumptions justified? What sets the 
seed field parameters? And how far can one trust 
the simulations themselves? 

Using vertically unstratified simulations performed 
by Sano et al. (2004), Pessah et al. (2007) provide 



a remarkably simple fitting formula for a, which 
states that its saturation value scales as the back- 
ground Bz/^fP^ where Bz is the mean strength 
of an assumed vertical magnetic field, P is the 
time-averaged gas pressure, and the simulation box 
length L is set equal to the pressure scale height h. 
Values of a - 0.01 correspond to B^J^ ^ 0.008, 
or equivalently a plasma /3 = SnP/Bz ~ 4 x 10^. Is 
such a background state realistic? Can one deter- 
mine analogous requirements based on simulations 
that assume a net toroidal field, and ultimately on 
simulations that include vertical gravity? 

4. The value of Am* . Is Am* truly and exactly 
100? What is the origin of this dimcnsionless num- 
ber? How sharp is criterion ([2]) in predicting the 
non-linear outcome of the MRI? According to our 
present calculations, criterion ^ for Am* ~ 10^ is 
only marginally satisfied at the rim and interior to 
the rim. We need to know the transport parameter 
a as a function of Am (see also point 3 above) to 
calculate more precisely the rate of gas flow from 
the rim. 

5. The value of Xmct, tot- How can one measure the 
abundance of gas-phase atomic metals in transi- 
tional disks? Our finding that the MRI is viable 
depends on having more than 1 out of every ~300 
metal atoms such as magnesium remaining in the 
atomic gas phase. These free metals provide the 
ions that coUisionally couple the neutrals to the 
MRI. 

6. Aerosols. What is the abundance of macro- 
molecules/grains having sizes ^lOA to ~0.1/im 
and is it large enough to interfere with the standard 
electron-driven chemistry we have presented? On 
the optimistic side, transitional disks are observed 
to have much less micron-sized dust than conven- 
tional T Tauri disks. We have calculated that frac- 
tional electron densities in the rim are large enough 
to withstand, to some degree, the effects of electron 
adsorption and ion recombination onto grains. On 
the other hand, radiation pressure cannot purge 
gas of small aerosols because of their short aerody- 
namic stopping times and low efficiencies for inter- 
cepting radiation. 

7. Accretion rate vs. X-ray luminosity. Any mature 
theory of X-ray-driven MRI should predict M as a 
function of Lx. But without knowing a (Am) (see 
points 3 and 4 above), we cannot derive this re- 
lation with confidence. At the very least M must 
increase with Lx- Is there observational evidence 
for such a trend? In fact, M is known to in- 
crease with Af* (e.g., Hartmann et al. 2006). And 
Lx is observed to increase with (Guedel et 
al. 2007; Telleschi et al. 2007). But frustratingly 

enough, when M is plotted against Lx, whatever 
trend might be present is obscured by large scatter 
(Telleschi et al. 2007). 

A shorter-term project would be to measure the 
hard. Ex ^ 3 keV spectra of individual sources. 
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The penetration depth of X-rays scales strongly 
with photon energy, as E^^^. 

8. From conventional T Tauri disks to transitional 
disks. We have studied the equilibrium dynamics 
of transitional disks but have said little regarding 
origin. Presumably transitional disks evolve from 
conventional T Tauri disks. Are transitional sys- 
tems, as a group, older? The age of TW Hyd is 
estimated to be approximately 10 Myr, which in- 
deed is old for a T Tauri star. 

How does the evolution from conventional T Tauri 
to transition system take place? Unfortunately, we 
cannot answer this question by merely extrapolat- 
ing our theory to earlier times. Our theory is de- 
signed to explain systems with large, well-defined 
holes, and would break down as Orim approaches 
the stellar radius, if only because rim gas would 
be ionized thermally rather than by X-rays. Data 
for conventional T Tauri disks clearly do not fit 
the relations plotted in Figure 4; their rim radii 



are on the order of ~0.1 AU but their accretion 
rates range up to ~lO~^M0yr~^. The origin of 
transitional disks will likely remain a mystery un- 
til we solve the decades-old puzzle of how conven- 
tional T Tauri disks, objects without large holes, 
accrete (see, e.g., Hartmann et al. 2006). The so- 
lution might involve some combination of the MRI 
plus gravitational instability (e.g., G96). 
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